Fractional revivals of the quantum state in a tight-binding chain 
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We investigate the time evolution of a Gaussian wave packet (GWP) in a tight-binding chain with 
a uniform nearest neighbor (NN) hopping integral. Analytical investigations and numerical simu- 
lations show that fractional revivals of the quantum state occur in this system, i.e., at appropriate 
times, a GWP can evolve into many copies of the initial state at different positions. The application 
of this quantum phenomenon to quantum information transfer in solid-state systems is discussed. 
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I. INTRODUCTION 

Fractional revivals of a quantum state occur when 
the wave function evolves in time to a state, at a spe- 
cific instant between two full revivals, which can be de- 
scribed by a superposition of states with equal ampli- 
tudes, each of which has the same shape as that of the ini- 
tial wave packet, yet with different spatial distributions. 
The fractional revival is an interesting phenomenon in 
quantum mechanics, which has no analog in classical 
physics. It has been studied extensively, both theoret- 
ically [E[2,i,i,[| and experimentally gSllli- Most 
of the studies have been devoted to continuous systems, 
such as the Coulomb potential and the infinite square 
well. In this paper, we will focus on the phenomenon of 
fractional revival in a discrete solid-state system. As an 
illustration, a simple tight-binding chain with a uniform 
nearest neighbor (NN) hopping integral is investigated 
analytically and numerically. Such a model is used to 
describe the Bloch electronic system in condensed mat- 
ter physics and now the qubit array relevant to quantum 
information applications due to its equivalence with the 
XY spin chain. It is found that fractional revivals of a 
Gaussian wave packet (GWP) occur in a discrete system. 
This opens up the possibility of performing high-fidelity 
quantum information transfer (QIT) and creating long- 
range entanglement by employing this novel feature of 
the time evolution in this solid-state system. 

Quantum information processing in solid-state systems 
has attracted widespread attention because of the poten- 
tial scalability of devices. Within this context, quantum 
state transfer (QST) from one place to another in such 
a system becomes a crucial issue and has been analyzed 
theoretically [13, M M M M M M MM,MMM- 
A great advantage of this approach is that no dynamical 
controls are needed after one prepares the quantum state. 
In one of the pioneering works [l^ , Bose considered a reg- 
ular one-dimensional spin chain with Heisenberg interac- 
tions, which is able to transfer a quantum state over a 
reasonable distance with the aid of a distillation process. 
Since then, a number of interesting proposals have been 
made for quantum communication through spin systems 
to improve the fidelity of the QST. One of them is to 
choose the proper modulation of the coupling strengths 



as suggested in [l6|, |l7|, [l^ . In such a system, although 
an arbitrary local quantum state will spread as the time 
evolves, after a period of time the dispersed amplitudes 
will "refocus" at the receiving location of the chain. So 
perfect state transfer can be realized. Another approach 
makes use of gapped systems. The advantage of these 
schemes is that the intermediate spins are only virtu- 
ally excited. In this case, the two separated qubits are 
coupled and realize the entanglement of two points pT| . 
This ensures that the transfer of the single-qubit state is 
achieved with a very high fidelity. 

It should be pointed that QST and QIT are two differ- 
ent concepts. QST means a local quantum state changes 
its location. It usually corresponds to the translation or 
refiection of the initial wave function. On the other hand, 
typical quantum information means the way (or mode) of 
the superposition of two orthogonal states. A complete 
QST can achieve perfect QIT. However, it is sufficient 
but not necessary. 

So far, almost all the proposed schemes for QIT are 
based on the fact that a properly designed qubit array 
can provide a unitary evolution operator U{t) which can 
accomplish the task of transferring the local state {tpA) 
at position A to the target state |V'b) located at B via 
the process of time evolution \%Pb) = U{t) \tpA)- 

The phenomenon of quantum revival is an example 
of this process, which transmits the complete local state 
over the distance if {tp b\U {t) a) ~ 1 (see Fig. [T^). 
In this case, we can say that perfect QST and QIT are 
both achieved. However, theoretically, if the data bus 
guarantees the partial revival, i.e., U{t) \i1^a) — ub + 
uc |'0c) with {ipB \tpc) = (see Fig. [TJd), the quantum 
information encoded in the target state {iPb) can be also 
extracted from the final state U{t) \4>a)- In this case, 
the quantum state can not be transferred perfectly, but 
perfect QIT is accomplished. We will show that such a 
scheme can be realized based on the fractional revival of 
a quantum state. 

In general, there are two ways to employ the tight- 
binding model as a data bus for the QIT: (1) The qubit 
array is usually described by a Heisenberg spin chain sys- 
tem [lo| . Within the context of quantum state transfer, 
only the dynamics of a single magnon are relevant. Thus 
in the single-magnon invariant subspace, this model can 
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hopping integral — J. The Hamiltonian can be written as 
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FIG. 1: (Color online) Schematic illustrations of quantum 
information transfer implemented by the revivals of a polar- 
ized electronic wave packet: (a) at instant t — Trev, the ini- 
tial wave packet at A is in full revival at B; (b) at instant 
t = Trevfi, it is in half revival at B. 



be mapped into the single spinless fermion tight-binding 
model. The quantum information is encoded in the su- 
perposition of the single and zero particle states. (2) On 
the other hand, the quantum information can also be en- 
coded in the polarization of the Bloch electron (Fig. [T|). If 
the spin state of the Bloch electron is a conserved quan- 
tity for the Hamiltonian of the medium, the spin state 
cannot be influenced during the propagation, no matter 
how the spatial shape of the wave function changes [23| . 
In this case, the locality of the final state is crucial. We 
will discuss the two schemes in detail in section IV based 
on the formalism of the fractional revival in the tight- 
binding model. 

This paper is organized as follows. In Sec. II the model 
is presented. In Sec. Ill we formulate the general theory 
for the fractional revivals of a GWP in the tight-binding 
model. In Sec. IV we show numerical results that sub- 
stantiate the analytical results. Sec. V is devoted to the 
application of the theory for QST and QIT. Sec. VI is 
the summary and discussion. 



II. MODEL HAMILTONIAN 

For completeness, we consider a non-interacting spin- 
1/2 fermion system on a one-dimensional lattice, which 
is a simple tight-binding TV-site chain with uniform NN 



(1) 



where ^ denotes the fermion creation operator at j-th 



site with spin a =tii- Because this Hamiltonian does 
not contain a spin-spin interaction, the polarization of 
the spin of an electron is not changed as time evolves. 
Then the problem about the transfer of the spin state 
can be reduced to the issue of charge transfer. In this 
sense, we can concentrate on the spinless fermion model 



H 



N-l 

i=i 



j+i 



h.c. 



(2) 



where Oj denotes the spinless fermion creation operator 
at j-th site, and the open boundary condition is applied. 
In section V, the spin degree of freedom will be reconsid- 
ered for the discussion of QIT. Introducing the Fourier- 
transformed fermion operator 



N 



where k = mr/ {N -I- 1), n = 1, 2 
(HI) can be diagonalized as 



-^sin(fcj)a], (3) 
N, the Hamiltonian 



H 



fc 

— 2 Jcos fc. 



(4) 



with the single-particle eigenstate 

N 



fc) = y27(ivTT)^sin(fcj)|j), 



(5) 



where 



^)=4lO), |j)=a]|0). 
In the large- TV limit, a discrete coordinate system ap- 
proaches a continuous one. Correspondingly, state ([5]) is 
a standing wave which is the eigen wave function of the 
infinite square well. What is more, in lower energy re- 
gion fc ~ 0, the spectrum is « — 2J(1 — fc^) ^ 2Jk^, 
which is very close to that of the infinite square well. 
This indicates that for low energies, the physics of the 
tight-binding chain is approximately the same as that of 
the infinite square well. Accordingly, for an initial state 
which can be expanded by the lower eigenstates, its time 
evolution should be similar to that of the infinite square 
well approximately, which is the typical paradigm to il- 
lustrate fractional revivals It implies that, for a wave 
packet with low energy, the well-defined fractional revival 
formalism should be valid in such a discrete coordinate 
system. 
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III. FRACTIONAL REVIVALS OF GWP 

In this section, we investigate the time evolution of a 
GWP in the tight-binding system analytically. Although 
the formalism presented in this paper has been well es- 
tablished for continuous systems, we derive it exactly to 
show that this formalism can be extended to discrete sys- 
tems for some special states such as GWPs. 



A. General formalism 

We consider an initial state in the low-energy range, 
which can be expanded by the lower eigenstates. A typ- 
ical state meeting this condition is the zero-momentum 
GWP, which can be expressed as 



1^ (^o)> = 



1 



N 



^a^ij~Nof/2 



(6) 



where fli = X^^i ^ " is the normalization fac- 

tor; No denotes the center of the GWP and 1 + A/2 < 
No < N - A/2 ensures that the entirety of the GWP is 
situated within the chain approximately. The factor a 
determines the half-width 



(7) 



of the GWP in real space and also the range of the spec- 
trum related to its eigenstate expansion. 

Firstly, we study the revival of a GWP in the frame- 
work of the spectrum-parity matching condition (SPMC) 
[13, HJ, which is the basis for the investigation of the 
fractional revival. Actually, using the Fourier transfor- 
mation ([3]), the GWP ^ can be written as 



IV' = 



1 



' sin(fciVo)e 



(8) 



where = '^k^^^^^i^^o)s~'' is the normalization 
factor. For a system which has mirror symmetry, we 
have [P, H] ~ 0, where P is the reflection operator de- 
fined as P\j) ~ I + 1 — j) . For the model Hamiltonian 

which has the mirror symmetry, the eigenstate 
k) = Pk 



k) with Pk = ±1. So the mirror 



satisfies P 

counterpart of the GWP ^ has the form 

P|V(iVo)) = |i^(A + l-Ao)) 



(9) 



1 



^(-l)"+^sin(fcA^o)e^ 



with the relation k = mr/ {N + 1). For the GWP ^ 
with a large enough a, it can be expanded by the eigen- 
state k ) with the eigenvalue Sk and parities pk in the 



following way: 



Ek = n^AE, 



Pk 



-1)" 



(10) 



where AE = 2 Jtt^/ (A'' -I- 1)^ is the greatest common di- 
visor of all the possible level differences. Obviously, the 
dispersion relation and the corresponding parity (fTO|l sat- 
isfy the SPMC, which leads to the following conclusion: 
for an initial state |(/)(Ao,0)) — \tlj{No)), at the instant 
t = Trev = 7r/A_B, it evolves into 



cI^{No,^))=P\cI){No,0)) 



(11) 



Here we define Trev as the revival time, at which the 
wave packet revives as the mirror refiection of the initial 
one. This is the same as in an infinite square well, and 
is usually called a full revival. 

In the following, we will show that fractional revival of 
the GWP also occurs in the tight-binding model. We 
start our investigation with the time evolution of the 
GWP \il^{No)) in the system At the time t, the 

initial state |0(Ao,O)) has evolved into 



\^iNo,t)) 



„i2Jt 



e ""|0(Ao,O)) 



(12) 



/a 



^ sin(fcAo)e 



-k^ /2a^-i2Jk^t 



Here, we take H ~ ^2J^j,(l — fc^)a|afc as the small- 
k approximation for GWP. Neglecting the overall phase 
gj2 jt have 



\<t>{No,t)) 



: ^ sin(A:Ao)e 



-k^ /2a^-i2Jk^t 



(13) 

The shape of \(j){No,t)) in real space as the result of 

the interference of the standing waves fc^ depends on 

the time t and A^o- In general, the feature of revival 
is characterized by the autocorrelation function A{t) = 
{(f>{No,0)\ (j){No,t)). In this paper, we choose another 
quantity, the mirror fidelity 



F(t) = ((/)(Apo,O)|0(Ao,i)> 



(14) 



where Npo — N + 1 — No denotes the mirror position of 
No , since the revival of the wave packet at different loca- 
tions is desirable for the task of quantum state transfer 
and the generation of entanglement. The fidelity is a 
quantity to characterize the QST. 

Now we focus on the special instants r = pTrev/q, 
where p, q are two mutually prime integers. In the frame- 
work of quantum information, F[t) — 1 indicates perfect 
quantum state transmission. At the moment r, the fi- 
delity of QST can be expressed explicitly 



\F{r)\ - 



E(-i)^ 



(15) 



where a„ = yJl/Vl2 sin(fcAo)e~'° is the expansion 
coefficient, which is the starting point of our discussion. 
Obviously, when p = we have r = T^-ev and F{t) = 1. 
This result is in agreement with the prediction from the 
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SPMC, and can be employed to perform the perfect QIT. 
As discussed in the introduction, it may be possible to 
accomplish QIT in the case that the quantum state is not 
transferred completely. 

In the following, we will demonstrate it in the frame- 
work of the well-defined formalism of fractional revivals. 
It is easy to see that exp(— zpn^Tr/g) is a periodic function 
with period I, i.e. 



g-ip(n+i)-7r/g _ g-ipn tt/ij 



where I is determined by q 
I = 



2q (odd q) 
q (even q) 



Performing the Fourier transformation, we have 

i-i 



where 



-ipn 7T / q 



l-l 



r=0 



^ — i2mzr /I 



(16) 



(17) 



(18) 



(19) 



n=0 



A straightforward calculation shows that br satisfies the 
relation 

K^bi^r, (r = l,2,...,Z/2-l). (20) 
Then at the instant r, the state \(t> {Nq, 0)) evolves into 



\^iNo,T)) 



^ ^ n r— 

bo\^P{No))~bi/^P\yj{No)) (21) 

l/2-l 

+ br[yj{N+)) + \^{N-))], 

r=l 



where = Nq±2{N +1) r/l. It is clear that, at the 
time T = pTrev/q, the initial GWP at Nq has evolved 
into / sub-GWPs. Each of them closely reproduces the 
shape of the initial one but at the positions Np^, 
(r = 0, 1, 2, 1/2 — 1) and the probability of the sub- 

GWP at 7V± (Npo) is \br\^ {\bi/2\^). Eq. ^ shows 
that the evolved state at time r is the superposition of 
as many as I "clones" of the initial wave packet with 
amplitude 6,- {r — 0, 1, • • • , Z/2 — 1) and different spa- 
cial positions. The final state \(f> {Nq, t)) is formed in the 
following process. First, split the original wave packet 
into I copies of the initial wave packet, each with the 
probability \br\^ . Then translate the copies to the posi- 
tions N^ = No±2{N + l)s/l, {s = 0, 1,...,V2). For 
the sub-GWPs with beyond the chain, i.e., N+ > N 
or N- < 0, we have \tp {N+)) = e'T |-0 {N+ - N)) and 
\i}{N-)) = e'"" \^p{\N-\)) from ©. This indicates that 



(N^)) will be reflected with 7r-shift when iV^ is be- 
yond the chain, and there are always two sub-GWPs at 
the initial position A^o and its mirror place Npo respec- 
tively. 

Defining the function 



fiNA,NB) ^ {cbiNA,0)\ cbiNB,0)) . 
thus the corresponding fidelity is 

F{t) = {<j>{Npo,0)\^{No,T)) 

l/2-l 



(22) 



(23) 



7/2 



J2 brfiNpo,N^). 



r=l,X=± 



Here, term {(j){No,0) \ (j){Npo,0)) has been ignored since 
we are only interested in the case with \NpQ — No\ ^ A. 
Furthermore, from the relation 



-i27T{r/l+p/q)^ 



where r' = r + 2pl/q, we have 



(24) 



(25) 



i.e., the probability of each sub-packet is l/q. However, 
Eq. (PT|) shows that the final I sub-GWPs may be not 
orthogonal due to the reflection and then their super- 
position determines the shape of the final state. On the 
other hand, the overlap of two neighbor GWPs should af- 
fect the shape of the final state if the sub-GWPs are too 
numerous. We will discuss this in Sec. V with the aid 
of numerical simulation. Actually, the above conclusion 
is also valid for an arbitrary initial state which satisfies 
the low-energy condition. 



B. Application of the formalism 

In the following, our aim is to apply the formalism 
developed above to several special cases, which should 
be useful for the transmission of quantum information. 
From the above analysis, we have |i^(T)| = 1 at r = T^e^j 
which shows that the initial GWP can be revived com- 
pletely at this instant. However, according to the for- 
malism, Trev is not the shortest period to perform per- 
fect state transfer. If a proper Nq is chosen, a shorter 
period can be obtained. Here we investigate the case 
with A^o = Af/3 to illustrate this point. When we con- 
sider the case p/ q = 1/3, the corresponding period of the 
Fourier transformation in Eq. (jl7p is Z = 6. Then at time 
T = T,.ei,/3, we have 



N T 
¥' ~3~ 



P 



(26) 



It shows that at r = T^evl'^ the initial GWP recurs to- 
tally at its mirror counterpart. 
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This result can also be explained in the framework of 
the SPMC. Actually, for the eigenstate expansion of the 
initial GWP \tp {N/3)), it is easy to find that the expan- 
sion coefficients of the levels n — 3m (m is the integer) 
vanish. Then the greatest common divisor of the effective 
levels for the state \tp (N/S)) is AEeff ^ 6Jn^/ {N + 1)^ 
= 3Aii^, and the corresponding recurrence period T/g^, — 
Trev/-^- A similar analysis can be applied to the case of 
the initial state \4'{N/2)). In this case, the expansion 
coefficients of the levels n — 2m vanish, and the corre- 
sponding recurrence period T^g„ = Trev/8. 

Except these special cases, the greatest common divi- 
sor for the effective levels of all GWPs at other positions 
remains AE. Nevertheless, if the number of the levels, 
which determines the greatest common divisor to be AE 
is few, an "effective" greatest common divisor should gov- 
ern the recurrence time dominantly. For example, if we 
take iVo — N/ni, (m > 3), we also notice large partial 
revivals at r = Trev/m. For Nq = 7V/4, at r = T^ev/^ 
the wave packet evolved into 



N T 



= (foo-fei) 



{bi-h2)P 



f ^) ) (27) 



which represents four cloned small GWPs at positions 
N/A and 3iV/4. Note that the superposition leads to |6o — 
= (2 - \/2) /4 w 0.146 and \hi - 62!^ = (2 -I- V2) /4 
« 0.854, which indicates that at the time t = Trev/^ the 
initial GWP splits into two sub-GWPs, and the one at 
NpQ is the large partial revival. 

Another example to illustrate this mechanism is when 
the initial state is a superposition of two GWPs with 
— Na, Nb respectively, i.e. 



\<^{Na,Nb)) = 



1 



-^[\^iNA)) + HiNB))] 



= J2 ^ [sm{kNA) + sin(fciVs)] ( 



The levels with vanishing expansion coefficients are de- 
termined by 



n{NA + NB) = 2Nk 
or n{NB - Na) ^ {2k - 1)N. 



(29) 



One of the solutions is Nb = 2Na = 2iV/3. The vanish- 
ing levels are n — 2k and n = 3{2k — 1) with fc = 1, 2, .... 
Then the greatest common divisor is 2AAE. The wave 
packet will be revived at multiples of Trev/2A for this 
case, as shown in Fig. ^p. 



IV. NUMERICAL SIMULATIONS 

The analysis above is based on the assumption that 
the spectrum of the system is quadratic. However, this 



is only approximately true in the lower energy range. In 
order to demonstrate the fractional revival in the tight- 
binding system and to show how exact the approximation 
is, in this section we will exhibit numerical simulations 
for finite size systems. 

We start our investigation from the revival and frac- 
tional revival in the general case. We consider the time 
evolution of an initial GWP with A = 24 and iVo = 50 
in the system with N = 500. According to the formalism 
we have Trev = {N + 1)^ /(27rJ) « 4.0 x lO'^/J, which is 
taken as the unit of time t in the numerical results. The 
fidelity |F(i)|^ over the interval t G [0, 6Tret,] is plotted 
in Fig. [2^. It shows that the fidelity has peaks around 



the instants r = T^-ev, 5Trev, 5Tre 



which is in agree- 



ment with the formalism in Sec. II. Interestingly, there 
exist many regular small peaks between two neighbor big 
peaks. We present the small peaks in Fig. [2]3 to show the 
details of the small peaks. According to the general for- 
malism, for local GWP, at instants r/T^ev = 1, 1/2, 1/3, 
1/4, 1/5, the corresponding I i^(r) I should be equal to 
the values of r/Trev The plot in Fig. ^Bp is in agreement 
with the analytical results with high accuracy. In order 
to demonstrate the mechanism of the fractional revival 
more explicitly, we calculate the profile of the evolved 
state 



\UNo,t)\ = \{i\ (/.(7Vo,t))| 



(30) 



at r/Trev as a function of the position i and plot it in Fig. 
[21 From the analytical results, at instants r/Trev — 1/5, 
1/4, 1/3, 1/2, and 1, the cloned sub-GWPs have the 
probabilities of 1/5, 1/4, 1/3, 1/2, and 1, which result 
in the maximum values of the corresponding sub-GWPs 
to be 0.089, 0.099, 0.114, 0.140 and 0.198. The numeri- 
cal results, the number and the shapes of the sub-GWPs 
are in good accord with the theoretical prediction ap- 
proximately. Now we turn our numerical investigation 
to the special cases. Let us look at a wave packet ini- 
tially localized at Nq = N/3. The result introduced in 
Eq. (pS)) shows that the initial GWP recurs totally at 
its mirror part at time r = Trev/S. The numerical re- 
sult in Fig. 2^ shows that the theoretical analysis is in 
agreement with the result of numerical simulation rep- 
resented by \F{Trev/3)\^ = 1, which has also been well 
explained from the SPMC. On the other hand, we also 
demonstrate the evolution of |$ {N/3, 2N/3)) in Eq. ^ 
numerically. In Fig. |4)d, it shows that the first revival 
time is around Trev/'^A which is in agreement with the 
analytical result. It also indicates that the proper choice 
of initial wave packet can revive in a shorter time, which 
can be used to transfer long-range entangled GWPs in 
the discrete system. Numerical simulation is also per- 
formed in the same system but with Nq = N/4. The 
theoretical calculation shows that such an initial GWP 
should revive at NpQ = 3A^/4 after the time Trev/ 4: with 
a relatively higher fidelity 0.854. The numerical results 
presented in Fig. [5^ and b depict the characteristics of 
the time evolution of the initial wave packet \il!{N/A)) 
via the fidelity \F{t)\ and the profile of the wave func- 
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FIG. 2: (Color online) Plots of the square of fidelity, |-F(i)j , 
for the initial GWP with A = 24, A^o = 50 in the system 
with = 500. 2(b) is a part of 2(a) over one revival time. 
The dashed line indicates that the square of the fidelities at 
T/Trcv = 1/2, 1/3, 1/4, ... are approximately in a line. 



tion |(/)i(7V/4, Tre„/4)|. From the analytical result of Eq. 
(PTjl . the behavior in Fig. [5^ and b can be explained. 
Actually, at the instant rreu/4, the initial GWP splits 
into four sub-GWPs, with two of them being at iVpo and 
two at TVq. The final shape of the wave function should 
be two cloned initial wave packets with probabilities of 
0.854 and 0.146 respectively. These will result in the two 
maxima of the wave functions, 0.183 and 0.076 around 
the positions 3iV/4 and Af/4. 

Based on the numerical results presented in this sec- 
tion, we conclude that the fractional revival phenomena 
for local wave packets can be observed in the discrete 
system. 




FIG. 3: The illustrations of the fractional revival via the time 
evolution of a GWP with A — 24, A'^o = 50 in the system 
with iV = 500. (a) r = 0, (b) r = Trev/^, (c) r = r,.e„/4, 
(d) r = Trev/^, (e) r = Trev/2, (f) r = Tr<,v It shows that 
at r = Trevp/q, the GWP splits into several sub-GWPs at 
corresponding positions with the same shape as the initial 
one. 



Heisenberg spin chain system [T^ . Within the context of 
quantum state transfer, only the dynamics of the single 
magnon is relevant. Thus in the single magnon invariant 
subspace, this model can be mapped into a single spinless 
fermion tight-binding model. The quantum information 
is encoded in the superposition of the single and zero 
particle states. (2) On the other hand, the quantum in- 
formation can also be encoded in the polarization of the 
Bloch electron (Fig. [T|). If the spin state of the Bloch 
electron is a conserved quantity for the Hamiltonian of 
the medium, the spin state cannot be influenced during 
the propagation, no matter how the spatial shape of the 
wave function changes p3|. In this case, the locality of 
the final state is crucial. We will discuss the two schemes 
in detail in the following based on the formalism of the 
fractional revival in the tight-binding model. 



V. FRACTIONAL FIDELITY OF QIT 



A. Scheme A: qubit array 



In general, there are two ways to employ the tight- 
binding model as a data bus for quantum information 
transfer: (1) The qubit array is usually described by a 



It is well known that, by employing the Jordan- Wigner 
transformation [15i|, the one-dimensional tight-binding 
chain with NN hopping is equivalent to a simple XY 
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1/3 2/3 

t/T 




FIG. 4; (a) Plot of the square of fidelity \F{t)\ , same as Fig. 
2b, but for A'o ~ N/3 over one revival time, (b) Plot of the 
square of fidelity over the time [0, Trev/4:] for the initial state 
which is a superposition of two GWPs with A'o ~ N/3, 2N/3 
respectively. It shows that for such an initial state, the revival 
period is reduced. 



chain. Such a system is usually used to depict the 
physics of the qubit array. In this paper, thejaasis of 
the XY chain is in the form Ylj with n — 0,1, i.e., 



and 



In general, the transmis- 



sion of a qubit state from the location A to B can be 
regarded as the following process. The initial qubit state 



an operation Uab to realize 



) is prepared at A. If one can find 
A 



Uab = u I'^^ + e^^^'^v 0^ 



(31) 



where fAB is the known phase for a given system, we 
say that the qubit state is transferred from A to B per- 
fectly. Then the perfect QIT is accomplished. Bose [l^ 
proposed that the operation Uab can be achieved in the 
qubit array by the time evolution of the system based on 



the fact that the saturated ferromagnetic state JJ, 




0.2 



0.1 



0.0 



(b) . 




A, 


1 



N/4 



N/2 

i 



3N/4 



N 



FIG. 5: (a) Plot of the square of fidelity \F{t)\ , same as Fig. 
2b, but for A^o = N/i over one revival time, (b) The profile 
of the evolved GWP at the time r = Trev/ 't- 



is an eigenstate of the model. Furthermore, it is found 
that perfect state transfer can be implemented if the sys- 
tem meets the SPMC. Unfortunately, for an array to sat- 
isfy the SPMC, it requires modulation [M [Ulil, [13 of 
the couplings between qubits, which is difficult to pre- 
engineer in experiments. 

Nevertheless, we can consider the transferred state to 
be not a single-qubit state at a certain site but a superpo- 
sition of single-qubit states localized in a small range of 
the coordinate space. Together with the saturated ferro- 
magnetic state, the quantum information can be encoded 
in such a single magnon Gaussian wave packet at iVo 



No 



/TTl 



This state contains the same quantum information as 
that of the single qubit state lip a)- So if the GWP §^ 
appears completely at another place, perfect QIT is ac- 
complished. In the following, we will show that if the 
GWP ([5^ appears partially, perfect QIT can also be ac- 
complished. 
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1.0 



FIG. 6: (Color online) Plot of the square of the fractional 
fidelity \Ff{t)f for the GWP over one revival time with A = 
24 in the system with N = 500 for (a) iVo = iV/6, (b) No = 
A'^/10. It shows that |_F/(f)|^ is close to 1 at several specified 
instants. 



According to the formalism of the fractional revival, we 
note that at t = pT/q there always exists a cloned sub- 
GWP of the initial state at the mirror position, with the 



probability lb. 



'l/2\ 



From the point of view of quantum 



information, theoretically, the information of initial state 
encoded in the initial state by factors u and v has been 
transferred to its counterpart completely, although the 
fidelity of QST F{t) may be far from 1. 

In order to depict this fact, we introduce the fractional 
fidelity, which is expressed as 



\Ffir)\ 



1 



0l/2\ 



\Fir)\ 

l/2-l 

-bi/2+ brfiNpo,N^) 

r=l.A=± 



(33) 



and is unity if the retrieved sub-GWP is the exact clone 
of the initial state. Then the QIT can be transferred 
perfectly, even the QST is not completely. 



0.0 



□ 



a 



□ 
A 



o 



N=700 
N=600 
N=500 
N=400 
N=300 



12 



16 



20 



24 



FIG. 7: (Color online) The square of fractional fidelity 
\Ff(Trev /'2)f versus A for systems with various lengths A'^ G 
[300,700]. It shows that enlarging A can enhance the frac- 
tional fidelity. 



To demonstrate the new definition of the fidelity and 
verify how much the initial state and the sub-GWP are 
alike, numerical simulation is performed in a finite sys- 
tem. The square of the fractional fidelity \Ff{t)\^ for 
the propagation of a GWP with A = 24 from Nq = iV/6, 
N/10 to Npo = 5iV/6, 9iV/10 in the system with N = 500 
within [0, Trev] is plotted in Fig. [51 It is obvious that at 
many instants r = TrevP/l, the square of the fractional 
fidelity approximately equal to 1 , while | J^/ (r) | is still far 
from 1 for many possible TrevP/Q- We also notice that 
the initial position affects the results strongly. These 
phenomena can be explained as follows. 

As discussed in Sec. Ill (A), the final state should 
be the superposition of the I cloned sub-GWPs. Al- 
though there always exists a cloned sub-GWP at the 
position Npo, the overlap of these sub-GWPs affects 
|F/(r)|. If the I cloned sub-GWPs are well separated, 
i.e., f{Npo,N^) = 0, then we have |F/(r)| = 1 directly 
from the Eq. ([55)) . On the other hand, the orthogonality 
of the I cloned sub-GWPs is determined by their num- 
ber and positions, which depend on the factors q and Nq. 
Obviously the smaller I is, the more perfect the fractional 
revival is. However, in the case that the initial position 
A'o and I satisfy the condition 



{N + 1- 2Nq)1 = integer x 2 (iV -f 1) 



(34) 



the fractional revival should be still perfect. The Eq. 
([Ml guarantees that the final state consists of several well 
separated cloned GWPs. This is in agreement with the 
numerical simulations. Fig. [5^ shows that |-F/(t)| ~ 1 
at T = Tre«p/12 (p = 1,2,- •• ,11), while it occurs at 
T = T,e.p/10 (p = 1, 2, • • • , 9) in Fig. B^. 
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B. Scheme B: Flying qubit 

There is another scheme for QIT when the spin de- 
gree of freedom is considered. We can imagine an elec- 
tronic wave packet with spin polarization as an analog of 
a photon "flying qubit" , i.e. a polarized photon where 
the quantum information is encoded in its two polariza- 
tion states. We define the solid-state "flying qubit", at 
the location A in a quantum wire, as the superposition 
of two orthogonal Bloch electronic GWPs |t)y^ and 
where 

\a =T, i)A = ^Y. e-"^(^-^-)^/^sV |0> • (35) 
j 

Obviously, the two orthogonal states evolve indepen- 
dently. Then an arbitrary state 

l*(0))A-"IT)A + Hi>A (36) 

will evolve to 

\^{r))^um,t))+v\<P{Ut)), (37) 

where |0(cr,t)) = Ej |0> with /(j,t) being a 
spin-independent function. The quantum information 
encoded in the spin state of the initial state (P5|) is car- 
ried along by the electron and unaffected by the transfer. 
Therefore, the initial state will be transferred to another 
location if f{j, t) is known and is still a localized function. 
Of course, the simplest case is that f{j,t) is a mirror or 
translation of the initial GWP [l^. In fact, if f{j,t) is 
partially local at several places, the quantum informa- 
tion can be transported to multiple receivers. This fact 
indicates that such a system can be used for a "quantum 
fanout" , which was recently proposed by A.D. Greentree, 
S.J. Devitt, and L.C.L. HoUenberg f25l|. In our work, we 
only employ a simple open chain without any dynamic 
control. In this sense, it acts as a solid-state based split- 
ter, entangler [26| , and quantum fanout. 

C. Validity of the schemes 

In Sec. Ill, the analytical conclusion is only valid for 
lower energy GWPs. For an arbitrary GWP, the factor a 
determines the behavior of the final state as time evolves. 
As pointed above, for the GWP with narrow width at fc ~ 
in momentum space, the effective dispersion relation is 
quadratic approximately. Then the width A of GWP 
in real space should affect the fractional fidelity of the 
fractional revival. Numerical simulation was employed to 
investigate the relationship between the fractional fidelity 
Ff{T = Trev/"^) and A with the GWP transferring from 



Nq = 50 to Nqp = N — 50. The numerical results for 
N = 300, 400, 500, 600 and 700 are plotted in Fig. [3 
It shows that when A tends to 24, F/(t — Trev/"^) ^ 1 
for different sizes of system. Thus it indicates that when 
fractional revivals in the discrete system are employed 
for quantum information transmission, the width of the 
chosen GWP should be more than 24 times the lattice 
spacing. 



VI. SUMMARY AND DISCUSSION 

In summary, we have studied the phenomenon of frac- 
tional revivals in a discrete system by theoretical analysis 
and numerical simulations of the evolution of a GWP in 
a tight-binding model. It is found that, for a proper 
chosen initial state, its fractional revival states have the 
same formalism as that in the infinite square well. On 
the other hand, numerical simulations show that the for- 
mulas of the theoretical analysis are very accurate for the 
GWP. 

We also proposed the concept of the fractional fi- 
delity Ff{t) when the fractional revival phenomenon is 
exploited to achieve QST in the solid-state system. We 
showed that the fractional fidelity approximately equals 
to 1 at many instants. With an appropriately chosen 
width, the GWP with a polarized spin state can be re- 
garded as a flying qubit in the solid-state system to im- 
plement quantum information transmission. 

It is worthwhile to discuss the applicability of the 
scheme presented above. Experimentally, the tight- 
binding model can be realized by a quantum dot array, 
SQUID array, etc.. In a real system, quantum decoher- 
ence is the main obstacle to the experimental implemen- 
tation of quantum information. In our scheme for quan- 
tum state transfer, the quantum decohcrcncc time lim- 
its the scale of the quantum wire. For coupled quantum 
dots, experiments show that the coupling strength J ~ 10 
meV [271 . Therefore, for a iV-site chain, the revival pe- 
riod Trev = {N + 1) V( Jtt) w 1.6 X lO^^TV^ On the 
other hand, the decoherence time of a quantum dot is 
T ~ Ims [28|. For example, considering a quantum dot 
array with N = 500, we have Trev ~ 4 x lO^^ms, which 
is much smaller than r. Moreover, if we perform n times 
full revival within the decoherence time t, the maximal 
size is 2.5 x 10^/^/n. Therefore, for n ~ 10^, the size of 
the system should be limited to lO'^, which implies the 
applicability of the scheme in practice. 
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